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Recently certain non-supersymmetric solutions of type Ilb supergravity were constructed [1], 
which are everywhere smooth, have no horizons and are thought to describe certain non-BPS mi¬ 
crostates of the D1-D5 system. We demonstrate that these solutions are all classically unstable. 

The instability is a generic feature of horizonless geometries with an ergoregion. We consider the 
endpoint of this instability and argue that the solutions decay to supersymmetric configurations. 

We also comment on the implications of the ergoregion instability for Mathur’s ‘fuzzball’ proposal. 


I. INTRODUCTION 

String theory has made great progress in understanding the microphysics of black holes. In particular, for certain 
(nearly) supersymmetric black holes, one is able to show that the Bekenstein-Hawking entropy S'bh = ^hor/dG, as 
computed in the strongly-coupled supergravity description, can be reproduced in a weakly-coupled D-brane description 
as the degeneracy of the relevant microstates [2] — for reviews, see [3]. The AdS/CFT correspondence [4] provides 
further insights into these issues by providing a dictionary relating the geometric description of the physics in the near¬ 
horizon region with the physics of a dual conformal field theory — see [5] for a review. In particular, the AdS/CFT 
indicates that Hawking evaporation should be a unitary process, in keeping with the basic tenets of quantum theory. 
The discussion of black holes in the context of the AdS/CFT correspondence makes evident that the path integral 
over geometries in the bulk may include multiple saddle-points, i.e., several classical supergravity solutions, as found 
e.g., in [6-8].^ Another point that was realized early on is that the geometric description of individual microstates 
would not have a horizon [9, 10]. 

In recent years, Mathur and collaborators have incorporated these ideas in a radical revision of the stringy description 
of black holes — for a review, see [11]. They argue that each of the CFT microstates corresponds to a separate 
spacetime geometry with no horizon. The black hole is dual to an ensemble of such microstates and so the black 
hole geometry only emerges in a coarse-grained description which ‘averages’ over the microstate geometries. In 
particular, this averaging should produce an effective horizon at a radius where the individual microstate geometries 
start to ‘differ appreciably’ from one another [12, 13]. Therefore in this scenario, quantum gravity effects are not 
confined close to the black hole singularity, rather the entire interior of the black hole is ‘filled’ by fluctuating geometries 
— hence this picture is often referred to as the ‘fuzzball’ description of black holes. The first support for this proposal 
came from finding agreement between the propagation time of excitations in the throat of certain microstate geometries 
and in the dual brane description [13, 14]. A further remarkable feature, that has drawn attention to these ideas, is 
that there is growing evidence that the microstate geometries may be smooth, as well as horizon-free.^ In the case of 
the D1-D5 system, smooth asymptotically flat geometries can be constructed corresponding to all of the RR ground 
states in the dual CFT [13, 15]. Despite their large degeneracy, this two-charge system will not produce a macroscopic 
black hole horizon. However, a large horizon can be produced by introducing a third charge, Kaluza-Klein momentum 
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[16, 17]. Recently progress has been made in constructing smooth microstate geometries in the D1-D5-P system [18- 
20]. While large families of such solitons are now known, a complete understanding of the three-charge case remains 
to be found. Further preliminary work on the four charge system of D1-D5-P-KK has also appeared [21]. 

In general, the preceding discussion connecting microstates with smooth geometries focuses on supersymmetric con¬ 
figurations. This raises the interesting question of how the fuzzball proposal would be extended to non-supersymmetric 
black holes. In particular, are there non-supersymmetric versions of the smooth horizon-free geometries corresponding 
to non-BPS microstates? Remarkably, Jejjala, Madden, Ross and Titchener [1] recently extended the known set of 
D1-D5 microstate geometries with a family of non-supersymmetric solutions, hereafter referred to as JMaRT solitons. 
The JMaRT solutions comprise a five-parameter family of non-supersymmetric smooth geometries which are asymp¬ 
totically flat.^ These solutions may be parameterized by the Dl-brane and D5-brane charges, the (asymptotic) radius 
of the internal circle with Kaluza-Klein momentum, and by two integers m and n which fix the remaining physical 
parameters. These integers also determine a spectral flow in the CFT which allows the underlying microstate to be 
identified. For m = n + 1, the JMaRT solitons reduce to supersymmetric solutions found previously in [15, 18, 19]. 

An important feature which distinguishes the JMaRT solitons from any of the analogous supersymmetric solutions 
is the presence of an ergoregion. As a consequence, in these non-supersymmetric geometries, there is an inner region 
(that extends to the origin) where states of negative energy are allowed. This then leads naturally to the question of 
whether or not the ergoregion produces an instability of the background. One possibility is that the ergoregion may 
lead to superradiant scattering which can produce a catastrophic instability in some situations [22-24]. However in 
the present case, this possibility is easily dismissed [1] because the solutions are horizon-free. Since the seminal work 
of Zel’dovich [25] on superradiant amplification of electromagnetic waves incident upon an absorbing cylinder, it has 
been known that the key ingredients for superradiance is the existence of an ergoregion and an absorbing surface. For 
black holes, the horizon plays the latter role but certainly the JMaRT geometries lack such a surface. 

Quite interestingly, there is another class of instabilities, which we simply refer to as ‘ergoregion instabilities’, that 
generically afflict spacetime geometries with an ergoregion, but no horizon. These instabilities were first discovered 
by Friedman [26], who provided a very general discussion. Explicit computations of the instability were later made in 
[27, 28] for the case of rotating stars with an ergoregion. There the existence of this instability was explicitly verified 
for a free scalar field in the background of a rotating star. According to Friedman’s general arguments however, 
the instability should also exist for electromagnetic and gravitational waves. Since the JMaRT solutions [1] have an 
ergoregion but no horizon, one might suspect that a similar ergoregion instability would arise in these geometries. 
The present paper then explicitly verifies the presence of an ergoregion instability for the JMaRT backgrounds with 
a variety of techniques. Further we consider the endpoint of the resulting decay and argue that it should be a smooth 
supersymmetric solution. 

Our results have immediate consequences for the endpoint of tachyon decay discussed in [29]. There, Ross extended 
the discussion of [30] to D1-D5 black strings for which he identified tachyonic string modes in a particular winding 
sector. He argued that the condensation of these tachyons would transform the spacetime to a JMaRT soliton. In 
conjunction with the above results, we see that these solutions cannot be the final endpoint of these decays but rather 
they should end with a supersymmetric microstate geometry. Our analysis and the ergoregion instability may also 
have interesting implications for Mathur’s fuzzball proposal more generally. 

The remainder of our paper is organized as follows: Section H provides a brief exposition on Friedman’s analysis 
[26]. In Section HI, we briefly review some of the features of the JMaRT solutions and present the main equations 
used in the subsequent analysis, namely the radial and angular wave equations for a free massless scalar field, as well 
as some of their properties. In Section IV we compute the details of the instability using a WKB approach [27]. We 
show explicitly that the instability exists for a general non-supersymmetric geometry of [1], and that it disappears 
for supersymmetric objects, as expected. In Section V, we use an alternative method, that of matched asymptotic 
expansions, to investigate the instability and its properties. The methods of sections IV and V are complementary, 
i.e., their regime of validity is different. We then perform a numerical analysis of the wave equation in Section VI to 
complement the analytical calculations. We find that the results of both analytical analyses agree remarkably well 
with the numerical results. In section VH, after summarizing the main properties of the ergoregion instability, we 
discuss various related topics: the endpoint of this instability; its consequences for Ross’s tachyon condensation [29]; 
general implications for the fuzzball picture of black holes. 


By considering orbifolding, this family can extended by a third integer [1] but we will focus on the original five-parameter solutions. 
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II. ERGOREGION INSTABILITIES 

There are two classes of instabilities that are of potential interest for the JMaRT backgrounds [1] (or non- 
supersymmetric geometries in general), namely: the superradiant instability, and the ergoregion instability. In this 
section, we demonstrate why superradiance is not present in these geometries, as first noted in [1], and we introduce 
the general argument of [26] that suggests an ergoregion instability is present. In the following sections, we verify the 
presence of the ergoregion instability with a complete analytic and numerical analysis of its properties. 


A. Geometries with an ergoregion and horizon: Superradiance 


For a general (stationary asymptotically flat) black hole, the equations describing spin-s fields can always be written 


as 


—Y + V(oj, r)'^ = 0 
art 


( 1 ) 


where iv was introduced with a Fourier transform with respect to the asymptotic time coordinate: 'l'(t) = e *‘^*5'(a;). 
The radius r* is a convenient tortoise coordinate and in general one finds: 


J r* ~ r , E ~ as r ^ oo , 

y e’’* ~ (r — r+)“ , E ~ (w — <I>)^ as r ^ r+ , 


where a is a positive constant. The potential <I> can be a rotational potential (in the Kerr geometry $ = mfl, with m 
an azimuthal number, and the angular velocity at the horizon) or a chemical potential (in the Reissner-Nordstrom 
geometry, $ = qQ, where q is the charge of the field and Q the charge of the black hole). 

For a wave scattering in this geometry, Eq. (1) yields the following asymptotic behavior:"^ 


«'i 


r(r-r+)-'“(“-‘^) asr^r+, 
7^e‘“’'+ e-‘“’' asr^oo. 


( 3 ) 


These boundary conditions correspond to an incident wave of unit amplitude from +oo giving rise to a reflected 
wave of amplitude TZ going back to +oo and a transmitted wave of amplitude T at the horizon — the boundary 
condition introduces only ingoing waves at the horizon. Assuming a real potential (which is almost always the case) 
the complex conjugate of the solution satisfying the boundary conditions (3) will satisfy the complex-conjugate 
boundary conditions: 


r T*(r —r+)'“^“ as r ^ r+, 

\ 7^*e-‘“’' -b e‘“’' as r ^ oo 


Now, these two solutions are linearly independent, and the standard theory of ODE’s tells us that their Wronskian, 
W = is a constant (independent of r). If we evaluate the Wronskian near the horizon, we get 

W = —2i(a; — <I>)|Tp, and near infinity we Hnd W = 2iw(|7^P — I). Equating the two we get 

= ( 5 ) 

LU 

Now, in general \'R\^ is less than unity, as is to be expected. However, for w — $ < 0 we have that \R\^ > 1. Such 
a scattering process, where the reflected wave has actually been amplified, is known as superradiance. Of course the 
excess energy in the reflected wave must come from that of the black hole, which therefore decreases. 

Superradiant scattering can lead to an instability if, e.g., we have a reflecting wall surrounding the black hole that 
scatters the returning wave back toward the horizon. In such a situation, the wave will bounce back and forth, between 
the mirror and the black hole, amplifying itself each time. The total extracted energy grows exponentially until finally 
the radiation pressure destroys the mirror. This is Press and Teukolsky’s black hole bomb, Hrst proposed in [22], 
This instability can arise with an effective ‘mirror’ in a variety of situations: a scalar Held with mass /r > w in a Kerr 
background creates a potential that can cause flux to scatter back toward the horizon [31]; infinity in asymptotically 
AdS spaces also provides a natural wall [32] that leads, for certain conditions, to an instability; a wave propagating 
around rotating black branes or rotating black strings may similarly find itself trapped [33]. 


Implicitly, we consider a massless field here but the discussion is generalized to massive fields in a straightforward way. 
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B. Geometries with ergoregion but horizon-free: Ergoregion instability 

Suppose now there is no horizon in the background spacetime. What changes with respect to the former discussion 
is the boundary conditions: since there is no horizon and no absorption. For this case, the boundary condition (3) at 
the horizon is replaced by some kind of regularity condition at the origin. We suppose the radial coordinate r now 
ranges from zero to infinity and we impose the following boundary condition: 

4' ~ Af{r) , r ^ 0, (6) 

where /(r) is some well-behaved real function. This ansatz encompasses for instance typical regularity requirements 
where, e.g., one chooses /(r) ~ with (3 > 0. Repeating the above calculation, one gets \TZ\'^ = 1. Therefore the 
absence of a horizon, which precludes any absorption, prevents superradiance and hence the superradiant instability. 

Nevertheless, geometries with an ergoregion but without horizons are the arena of another class of instability. This 
ergoregion instability was discovered by Friedman [26]. Even though his discussion was made in four-dimensions only, 
it is trivial to extend it to any number of dimensions. The instability arises because of the following [26]: Given the 
test field energy-momentum tensor we can associate a canonical energy 

Ss=[rT,^dSt, (7) 

Js 

where is the background Killing vector which generates time translations in the asymptotic geometry. Now, 
because is space-like within an ergosphere, initial data can be chosen on a Cauchy surface S which makes Es 
negative. Moreover, it is shown in [26] that the energy can be negative only when the test field is time dependent. 
Then, since the field is time dependent, and since only positive energy can be radiated at future null infinity, the value 
of Es can only decrease further from one asymptotically null hypersurface S to another, say. S', in the future of S. 
Thus, the energy Es will typically grow negative without bound. This instability was computed analytically using a 
WKB approximation in [27] for rotating stars. There it was shown that the instability timescale is usually very large 
(typically larger than the age of the universe). The analysis of [27] was improved in [28] where further details of the 
instability were computed numerically. 

A key assumption above is that the system can not settle down to a negative energy configuration which, while time 
dependent, is nonradiative. Friedman [26] was able to rule out such marginal cases where Es is negative but constant 
for a four-dimensional massless scalar or electromagnetic fields. However, in fact, one is able to identify negative 
energy bound states for the JMaRT backgrounds — see Appendix D — and so a more thorough analysis is called for. 
Hence, in the following, we apply a variety of techniques to explicitly show that these microstate geometries suffer 
from an ergoregion instability. 


III. FORMALISM 


We now consider wave propagation of a free massless scalar field in the JMaRT backgrounds [1], and from this 
identify an ergoregion instability for these geometries in the subsequent sections. The JMaRT solutions are described 
in detail in [1] and are quite involved. We will provide a brief discussion of some of the properties of these solutions 
here, but will refer the reader to [1] for the full details. 

The JMaRT solitons are solutions of type Hb supergravity corresponding to three-charge microstate geometries of 
the D1-D5-P system. The system is compactified to five dimensions on T"^ x with the D5-branes wrapping the 
full internal space and the Dl-branes and KK-momentum on the distinguished S^. The notation is best understood 
by considering the construction of these solutions. One begins with the general solutions of [17] which contain eight 
parameters: a mass parameter, M; spin parameters in two orthogonal planes, ui, 02 ; three boost parameters, Ji, Js, Sp, 
which fix the Dl-brane, D5-brane and KK-momentum charges, respectively; the radius of the S^, R; the volume of 
the (which plays no role in the following). The geometry is described by the six-dimensional line element written 
in Eq. (2.12) of [1], which is parameterized by a time coordinate t; a radial coordinate r; three angular coordinates 
9, (j), tp; and the coordinate on the S^, y. 

One then imposes a series of constraints to ensure that the solutions are free of singularities, horizons and closed 
time-like curves. In particular, one focuses on a low-mass regime, < (oi — 02 )^, in which no black holes exist. 
Then one finds solitonic solutions where an appropriate circle shrinks to zero at the origin and the constraints ensure 
that this happens smoothly. First, M and R can be fixed in terms of the remaining parameters — see Eqs. (3.15) 
and (3.20) of [1]. Two quantization conditions constrain the remaining parameters in terms of two integers m,n [1]: 


j+j 


1 


s + 



= m — n , 


= m + n, 


( 8 ) 



5 


where j = . — < 1 and s = 


_ /£i_s^ ^ ^ -yy using the notation here that ci = cosh(5i and si = sinh(5i. Without 

ai — Y cicscp — o I I t I 

loss of generality, one assumes oi > 02 > 0 which implies m > n > 0. We also note here that the special case 
TO = n + 1 corresponds to supersymmetric solutions. 

This leaves a five-parameter family of smooth solitonic solutions. We can think of the independent parameters as 
the Dl-brane and D5-brane charges, Qi,Qd; the (asymptotic) radius of the y-circle, R; and the two integers, to and 
n, which fix the remaining physical parameters as [ 1 ] 


Qp = nm 


Q 1 Q 5 

i ?2 


J0 = -TO 


Q 1 Q 5 

R 


Jib — ^ 


QiQb> 

R 


(9) 


Of course, depending on the specific application, it may be more appropriate and/or simpler to describe the solutions 
using a different set of quantities. In our case, when we make explicit calculations of the ergoregion instability, we 
will fix the parameters n,TO, ai,ci and C 5 or Cp. As we are interested in non-supersymmetric backgrounds, we also 
impose to > n -I- 2. To conclude our discussion of notation, we add that the roots of r+ and r_, will also appear 
in the following but they are determined by M and the spin parameters — see Eq. (3.2) of [1]. 

The key ingredient producing the instability in the JMaRT solutions is the existence of an ergoregion. To verify 
the presence of the ergoregion, one takes as usual the norm of the Killing vector V = dt and using Eq. (2.12) of [1], 
calculates 




U 


f-Mcl 


( 10 ) 


where /(r) = + of sin^ 0 -|- a| cos^ 6 > 0 and Hi{r) = f{r) + Ms^, i = 1,5. It is then clear that V = dt becomes 

space-like for /(r) < M and thus an ergosphere appears at /(r) = M. An inspection of the metric also allows one to 
conclude the geometry rotates along (j), ip and y since gtrfi ^ 0, g^if, ^ 0 and gty ^ 0. The supersymmetric limit of the 
JMaRT solitons corresponds to take the limit M ^ 0 and 6i ^ 00 , while keeping the other parameters fixed, including 
the conserved charges Qi = MsiCi [1]. So, in the supersymmetric limit the norm becomes |Kp = —//\/ HiH^, which 
is always negative and thus the ergoregion is not present. 

Now consider the Klein-Gordon equation for a massless scalar field propagating in the JMaRT geometries, 

Implicitly, we are using the string-frame metric in which case one can think of eq. (11) as the linearized equation 
of motion for the Ramond-Ramond scalar. As described above, these backgrounds can be thought of as special 
cases of the general D1-D5-P solutions found earlier [17] and so one may apply separation of variables following [34]. 
Introducing the following ansatz ® 


d* = exp 


t , y ■ 

-IW-- - iX— + im.^ip J- im^ 
R ^ 


R 


Xi0)h{3 


( 12 ) 


one gets an angular equation 


1 


sin 29 dO 
and a radial equation® 


sin 29 


d9 


A- 




cos2 9 sin^ 9 


+ ^ (Qi 9 + al cos^ 6 *) | X = 0 , 


l_d 

r dr 


g{r) d 


r dr 


- Ah + 


(w^ — A^) , 2 ,.2 2\ r ^ ^2 

(r J- Msi J- Ms^) J- (oJCp J- Asp) 


R 2 


-irl-rl) 


(A — nTO .0 + 

(p 2 _ p 2 ^ 


+ (r+ - ) 


{ujg + Xd — nm^ + mmtjjY 


h = 0. 


(13) 


(14) 


^ Note that the negative sign for A corrects a typo found in [1] 

® Note the factor (r^ — t?.) that appears in the two last terms of the Ihs of (14), which are necessary for dimensional consistency, corrects 
the typo appearing in Eq. (6.4) of [1] 
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where g{r) = (r^ — r^)(r^ — and we used \/—g = r sin 6 cos 9\/H 1 H 5 (this is the determinant of the metric (2.12) 
of [1]). If we introduce a dimensionless variable 


X = 



(15) 


we can rewrite the radial equation in the form 


with 


and 


dx[x{x + l)d^h] 


•*21 


K^x + 1- v'^ + 


a; + 1 


/i = 0 , 


t .2 _ t -2 

.^2 _ ,2 \ 2 ', ' + ' - 

K - (u - A )-^2- ’ 

^ = ojg + A-d — m^n + m^m , 
C, = \ — m^pn + mp,m , 


Q = 


„2 „ 2„2 „2 „2 „2 

CiCsCp 


c|c|-s?s 2 


S1C155C5 
2 2 
l'^5 


S1C1S5C5 


^pC-p 5 


— 1 + A- — —(r^ + Ms^ + Msg) — {iocp + ^Sp)^ — 


ur — A 


M 


(16) 


(17) 


(18) 


The quantities w, A, are all dimensionless — the last three being integers. Again, we refer the reader to [1] 

for a detailed account of the quantities appearing above. The reader should take note that our notation is not in 
complete accord with that of [1]. That is, to simplify our formulae in the following, we have defined n = 1/cr, the 
inverse of the quantity a used there. 

Of critical importance in characterizing the solutions of the scalar wave equation is the sign of The term xi^ 
dominates at large x, determining the asymptotic behavior of the solution. In this paper we will mainly be interested 
in outgoing modes so we choose k} to be positive. The two remaining possibilities: = 0 and k} < 0, will be 

considered in the appendices. 

The angular equation (13) (plus regularity requirements) is a Sturm-Liouville problem. We can label the corre¬ 
sponding eigenvalues A with an index I, A(lo) = Aimiuj) and therefore the wavefunctions form a complete set over the 
integer /. In the general case, the problem at hand consists of two coupled second order differential equations: given 
some boundary conditions, one has to compute simultaneously both values of uj and A that satisfy these boundary 
conditions. However, for vanishing of we get the (five-dimensional) flat space result, A = I {I+ 2), and the associated 
angular functions are given by Jacobi polynomials. For non-zero, but small ^2 Qj we have 


A — l{l + 2.) + O ^ . 


(19) 


The integer I is constrained to be I > |m^| -I- \mp\. We will always assume af <C max(m^, m|) (with i = 1, 2), 
thus A ~ l{l + 2). Making this assumption implies we may neglect the terms proportional to Oi in the angular 
equation, but given the way A and u appear in the radial equation, the corrections to A may not be negligible when 
we determine lo. To ensure that fixing A = l{l + 2) is consistent in both the angular and radial equations we must 
additionally require 

< max (|r^-I-M(s^-I-Sg)!, MCp) , (20) 

so that the contribution to ly from the Oi dependent corrections of A are negligible (see (18)). 

Taking the complex conjugate of Eq. (13) we can see that the exact solution to the angular equation has the 
symmetry 


Aim{-uy*) = A*^{iv). 


(21) 
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With this symmetry, one can also check the following: 


{v^Y{u:,X)=v\-u:*-X), 

( 22 ) 

(^^)*(a;, A,m^,TO 0 ) = ^^(-w*,-A,-m^,-m^), 

(23) 

(C^)*(A, m,p,mY = C^(-A, -m^, -m^). 

(24) 


Therefore, from the wave equation (16) it follows that if w is an eigenvalue for given values of with 

eigenfunction h, then —to* is an eigenvalue for with eigenfunction h*. Furthermore, if is 

outgoing unstable, so is Since the symmetry simultaneously flips all the signs of to^,TO 0 ,A, without loss of 

generality, we can only fix the sign of one, e.g., TZe{uj) < 0. 

To conclude this section, we point out that the angular equation (13) can be recast in the somewhat more familiar 
form: 

1 ^ ( ■ o o , 

sm 6 »cos 6 »— + 

Sin 0 cos 0 dO \ dO J 

where 

A = A H- — — al. (26) 

This is just the equation for a five-dimensional scalar spheroidal harmonic [35] which arises, e.g., in the separation of 
Klein-Gordon equation in the background of a five-dimensional rotating black hole [36]. 


A- 




-A" 


cos^ 6 sin^ 9 


i?2 


/ 2 2 \ 2/1 
(a 2 — ai) cos 0 


X = 0, 


(25) 


IV. WKB ANALYSIS 


We now explicitly show that the JMaRT geometries [1] suffer from an ergoregion instability. As described above, 
this instability is due to the fact that the geometry has an ergoregion but no horizon. We shall identify modes of the 
scalar field that are regular at the origin, represent outgoing waves at infinity and grow with time. In this section, 
we follow the WKB analysis of [27] and show that it applies to the non-supersymmetric JMaRT solutions, with the 
same qualitative conclusions. 

To begin, we want to write the radial wave equation in the form of an effective Schrodinger equation. In order to 
do so, we first transform to a new ‘wavefunction’ H defined with 

Hx) =—r^==Hix). (27) 

\/a:(l-f x) 

Inserting this in (16), we get 

-d^H + Ues H = 0 , (28) 


where 


Ueff 


+ (1 — + k'^)x'^ + {I — v'^ + — C‘^)x -f I — 

4x2(1-Pa;)2 


(29) 


Now in order to simplify our analysis, we choose: A = 0, = 0, and large m^. With A 0, the waves see a 

constant potential at infinity and thus the amplitude of the outgoing waves can be suppressed there. We also consider 
I = m., 1 , modes, which are expected to be the most unstable. Modes with I ^ m.,p must be similar to modes with 
= 0 for some I and these are not unstable. With these choices, we have 


2 2 

2 2 ^+ ~ /-2 2 2 t -2 2 2,2 2,0 

K = oj —— , Q =n 4 = m + u: g + 2 wpmm^ , 


(30) 

(31) 


Instead of working directly with the frequency of the wave, it will be convenient to work with the pattern speed along 
the 'tp direction, which is the angular velocity at which surfaces of constant phase rotate. This velocity is proportional 
to 


U) 



(32) 



where the proportionality constant i? ^ is always positive. It is important to compare the sign of the pattern speed 
along ip with the sign of the angular velocity of the geometry along if) defined as usual by^ 




9tTp 


2Msp Cp cos^ 0 Rjn 
\/ H1H5 9 ipip 

2Qp cos^ 0 cos^ 6 R/n 
\/ H 1 H 5 9'i’i’ 


(33) 


where Qp = MspCp is the Kaluza-Klein momentum charge. So, when is negative, the wave is propagating in the 
same sense as the geometry. 

Now it is useful to introduce the polynomial 


V = Bx^ + (A + B)x^ + {q^ + A)x , 

which is positive definite in the range of interest (positive x). We also define 

^ C/eff ^_rl+M{sl + sl + ^p) ^_tI- 

m I ’ - i ?2 

Then, we can write the effective Schrodinger equation (28) as 

dlH + mlTH = 0, 


(34) 


(35) 


(36) 


with 


T = 


4x^(1 + x)2 


y2 


2 gmx 


'Bih — 


— x{w? — — 1) + ; 


(37) 


where we have dropped certain small contributions to T.® Now it is straightforward to factorize the potential T and 
write it in the form 


with 


T = 


4x^(1 


r(s^-y+)(s^-y_) 


y Orr^^ 

v 


/ gmx\ 
V V ) 


2 x^ — x{m? — — 1) + 


(38) 


(39) 


For general m , n the behavior of the potentials V+ and V- (see Fig. 1) is exactly the same as the one studied in 
[27], so we do expect an instability to arise, as will be shown below. However, and this is a key point, for the case 
m = n + 1 which is the supersymmetric case, we have 


gmx 
V+ = —^ 


gmx\^ 

} 



TO = n + 1, 


(40) 


which is always positive. Thus, this WKB analysis indicates that the supersymmetric solutions are stable, as expected. 

Hence our radial equation has been reduced to the Schrodinger form (36) with an interesting potential (38), which 
depends on the pattern speed (32). Now the problem becomes to tune this potential by adjusting in order that 
a ‘zero-energy’ solution can be found with the appropriate boundary conditions: regular at the origin and outgoing 
waves at infinity. Note that in a region where is above V+ or below V- (allowed regions), the solutions have an 
oscillatory behavior. In those intervals where is in between the curves of and V- (forbidden regions), the 
solutions have a real exponential behavior. 


^ Note that the geometry rotates simultaneously along the ip, <p and y directions. We find using of (2.1), (3.17) and (3.19) of [1]. 

® More precisely, we have dropped a term l/(mJ'P). This remains a very good approximation in the high-m,^ limit in which we are 
working. As an example, for n = 10 and = 10 the factor that we dropped is 10“^ smaller than the last term of (37). 
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V 



FIG. 1: Qualitative shape of the potentials V+ and V- for the case in which an instability is present. An example of data that 
yields this kind of potentials is (m = 14 , n = 10 , ai = 32 , ci = 5 , Cp = 5). The unstable modes are those whose pattern speed 
is negative and approach the minimum of V+ from above. Thus, they are nearly bound states of the potential well in V+ 
that can however tunnel out to infinity through V-. Choosing A = 0, the potentials V+ and V- approach zero as a; —> oo, which 
makes a tunnelling through V- easier. 


We proceed following [27] and study the scattering of waves in the effective potential constructed above. Consider 
a wave that comes from infinity with an amplitude Cin, scatters in the ergoregion and returns to infinity with an 
amplitude Cout- In particular, we introduce the scattering amplitude defined as 

(41) 

^in 

The presence of a pole in S {i.e., of a resonance) signals the existence of an instability. Indeed, a pole in S occurs 
when Cin = 0 and Cout 7 ^ 0, and this means that we have finite outgoing radiation for zero incoming radiation. Near 
the pole frequency Wp, the scattering amplitude can be written to lowest order as [27] 


5 - 


g*2i5o 


LU — UJp ' 


(42) 


where So is a constant scattering phase shift and tt;* is the complex conjugate of Up. Note that this expression 
guarantees that when the frequency of the wave is real, one has S'((x>)[S'(ti;)]* = 1 as required by energy conservation. 
Generically, we can write the pole or resonant frequency as 


UJp = OJr + i/t , 


(43) 


where Wr and l/r are, respectively, the real and imaginary parts of tOp. With this convention, a mode with positive r 
represents an instability, and r < 0 represents a damping mode, since the time dependence® of the resonant wave is 
given by We can then write 


UJ — LOr — i! T 


(44) 


To relate the amplitudes Cin and Cout we apply a WKB analysis. As we shall learn later on, the unstable modes are 
those whose pattern speed is negative and approaches the minimum of from above (see Fig. 1). The scattering 
problem has then four distinct regions, namely: I, the innermost forbidden region (0 < a; < xq); II, the allowed region 
where V+ is below (xq < x < xi); III, the potential barrier region where V+ is above E^ (xi < x < X 2 ); and finally 


® Our conventions differ slightly from those of [27]. There waves carry a time dependence while we follow [1] which introduces the 
separation ansatz (12) with a time dependence 
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the external allowed region where is below V- (x2 < x < oo). The unstable modes are those that have < 0 . 
Thus, they are nearly bound states of the potential well in that can however tunnel out to infinity through IT. 
In region I, the WKB wavefunction that vanishes at the origin x = 0 is 


Hi 



( 45 ) 


where Ci is an amplitude constant. Then, the usual WKB connection formulae and WKB wavefunctions allow us to 
relate Hi with the wavefunctions of the other regions and, in particular, with the incoming and outgoing contributions 
of the wavefunction Hiy in region IV, which can be written as 

r 

-im^ / 

J X2 

The WKB analysis yields the relation between the amplitudes Cq, Cy and Ci (see Appendix A): 


Vt da 


( 46 ) 


Hiy- 


Ce 




exp 


i 


Vt da 


Cr 


my"ri/4 


exp 


where 




1 

2 
1 
2 


(2r;+^)c'6 + z(2r;-l)c7 


7 = 



Vt dx —, 

4 


( 47 ) 


( 48 ) 


In ?7 = / VW\ dx ■ ( 49 ) 

The identification of the ingoing and outgoing contributions in ( 46 ) depends on the sign of Sy,. Indeed, one has 
4 ' oc e“*“*i 7 iv(a;)- If is negative the term represents the ingoing contribution, while the term 

describes the outgoing contribution (if > 0 , the terms proportional to and Cy in Hiy{x) 
represent, respectively, the outgoing and ingoing modes). Henceforth we consider the < 0 case (since this will be 
the unstable case), for which the scattering amplitude can be written as 


Cy _ 1(4772 - l)e*'^ + (4772 + l)e-*'^ 

Ce (4772 + l)e*'>' — 7(4772 — l)e“*4' 


The resonance peaks in the scattering amplitude occur at a frequency lon for which e + ie*'*' = 0 , i.e., when 
'Y(uj) = 7 jv where 


iNiuJN) = Ntt + - 


( 51 ) 


with N being an integer usually referred to as the ‘harmonic’. The easiest way to see that the resonance peaks must 
be near these (real) frequencies is to note that S'(77v) = —i while for 77 ^ 00, one has 5(7 ^ 7^7) = +i. So when 
77 ^ 00, one has generally £'(7) = +7, but when 7 = a peak occurs that changes the value of S from +7 to —i. 

We can now do a Taylor expansion of the functions that appear in S around 7 = jn. Defining 



iO—CON 


d 



Vt dx 




the scattering amplitude can be written as 


( 52 ) 


-a{ui - ujn) + 4 ^ - f 

a{uj - ujn) + 4 ^ 

-a{uj - ujn) + 4 ^ + f 

a{uj - ujn) + Vp- 


(53) 
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which, using (1 + i)/(l — i) = i, can be cast in the form 


S'- 




(54) 


This result takes the form (44). Hence the discrete spectrum of resonance frequencies wat is selected by condition 
(51). Further comparing (44) with (54), one has that the growth or damping timescale is given by 

T = Arj^a . (55) 


Now, a defined in (52) is always positive since as increases so does T and 7 defined in (48) (the area of the region 
in between the line and the V+ curve, and in between line and the V- curve both increase when Y^ increases). 
So, we are guaranteed to have a positive r and thus the negative Y^ modes are unstable. If we redo the computations 
to consider the Y^ > 0 case, the only difference is that in (46) the ingoing and outgoing waves are given instead by 
the terms proportional to C^ and Cq, respectively. This changes the scattering amplitude from S to S~^ and thus r 
to —T implying that the positive Y^p modes are damped. 

Though the resonance frequencies and growth timescales can be computed with numerical methods from (51) and 
(55), as we shall do in Section VIB, we can still make some further progress analytically by approximating the well 
of V+ by a parabola. Near the well, the potential V+ behaves generally as 


rn---1- On 

m 


(56) 


with Om < 0. The boundaries Xq and Xi are the roots of Yp, — V+, namely: Xq = x^ — [P^niYp, — and 

xi = x^ + [Pm{Ypj — Since VT vanishes at these boundaries one has 


a = 



dVr 

dY^ij 


dx. 


(57) 


Moreover, near the bottom of the well, only Yp, — V+ varies significantly with x, and we can assume that all the other 
quantities that appear in the integral of a are approximately constants given by their value at x = Xm (the accuracy 
of this assumption increases as Y^p approaches Om)- One then has 


V A I QTTtXrr 


^/ViXm) 


\jYpj — W(Xni) 2Xni(l “t“ Xni) 

with V+ given by (56), which yields for a the value 


[S,/, - V+] = dx , 


= tta/^ 


Y^ + 


Qmx„ 
P{Xm)\ 


[S^-y_(x,n)] 


- 1/2 \/PiXm) 

2Xni (1 “t“ Xni) 


(58) 


(59) 


Let us illustrate the use of the WKB method we have described in this section to compute the instability parameters 
in a particular configuration. Take, 


m = 14 ; n = 10 ; oi = 32 ; ci = 5 ; Cp = 5 ; (60) 

A = mp = 0; I = mp = 10 . (61) 

By approximating the well in by a parabola, as in (56), we get 

= -0.17894 ; x^ = 9.1537 ; P^ = 2759.4 . (62) 

The resonant frequencies are those that satisfy condition (51) with 7 (w) given by (48). For the fundamental harmonic 
{N = 0), we get 


= -0.173. 

The growth timescale of the instability is given by (55) with r?(wAr) given by (49). Again, for = 0 we get 


(63) 


T ~ 10^". 


(64) 



12 


Independently of the parameters of the geometry, we note that as grows, approaches Om, the value of the 
at its minimum. For the particular geometry parameters described in (60) we have (for A = = 0): 

= 10: S^ = —0.173, 

= 20: S^ = —0.176, 

= 40 : = -0.177. (65) 

This feature can be proved analytically, as was done in [27]. 

Let us verify consistency of our results. We have assumed that af ^2 <C 1 in order to do the approximation 
A ~ l{l + 2). Now, for the cases we listed above one has of ~ 10“^, which is inside the range of validity 

for our approximations. A different combination of parameters yields a different instability timescale, and resonant 
frequency, so there are geometries more unstable than others. The following refers to the fundamental harmonic, and 
are computed within the parabolic approximation. We work with the following parameters: 


Cl = 5 ; Cp = 5 ; A = = 0 ; I = = 10 . (66) 

For the fundamental harmonic we then get 

TO = 1400 ; n = 10 ; m = 32 ; E^ = -0.40502 , r ~ 3 x 10®^ (67) 

m = 12 ; n = 10 ; oi = 32 ; E^ = -0.104 , r ~ 3 x 10 ^^ (gg) 

m = 14 ; n = 10 ; ai = 3200 ; E^, = -0.1728 , t ~ 3 x 10^® (69) 

TO = 3 ; n = 1; ai = 32 ; E^ = -0.0148 , r ~ 3 x 1.7 x 10^^ (70) 


It also evident that the instability is much stronger for small values of to^, where the WKB is expected to break 
down. 

To conclude this section, let us consider the regime of validity of the WKB approximation with more detail. A 
standard analysis of Eq. (28) suggests the WKB approximation is valid for \dxUes\ ^ |C4ffp, which can be rewritten 
as \dxT /T'^{x)\ <C to^. So, for large to^, the WKB approximation seems to be valid quite generally. However, we 
must sound a note of caution. As we already remarked, Eq. (65) shows that as grows, E^ approaches am, the 
value of the V+ at its minimum — this can be proved analytically [27]. So when to^ becomes very large, the two 
turning points are very close and the WKB analysis breaks down because T{x) 0. So we conclude that the WKB 
approximation used in this section should be valid in a regime with large m^, but not exceedingly large. In any event, 
it is clear that the instability is strongest for small values of to^, when the WKB analysis is certainly not valid. So, 
in the next two sections we will compute the features of the instability using complementary methods valid for small 
values of m^. We will also find remarkable agreement between all three approaches. 


V. MATCHED ASYMPTOTIC EXPANSION ANALYSIS 

The WKB analysis described in the last section appears to be strongest when describing solutions for which 
K~^ ~ C)?) but in general this corresponds to solutions with high angular momentum. In the sense that the timescale 
of the instability due to these modes is largest, they are the least unstable. Conversely, the matched asymptotic 
expansion that we use in this section becomes valid when k~^ > C.,^, they are the dominant decay modes. As an 
additional bonus, the eigenvalues are determined explicitly through algebraic constraints. Having both approximations 
at our disposal allows us to accurately calculate the eigenvalues for most of the allowed parameters. 

We follow a matching procedure introduced in [37], which has previously been used for studying scalar fields in 
three-charge geometries by Giusto, Mathur and Saxena [20], in the JMaRT backgrounds [1] and also in [24, 31-33, 38]. 
The space is divided into two parts: a near-region, a; <C /3, and a far-region, x ^ a, such that a (3- The radial 
equation is then solved approximately and the appropriate boundary conditions applied in each of the two regions. 
Finally, we match the near-region and the far-region solutions in the area for which they are both valid, a <C r <C /?. 
This gives a set of constraints, the solution of which gives the eigenvalues. Performing this analysis for the radial 
equation (16), we shall see that the only solutions which are regular at the origin and purely outgoing at infinity are 
finite as a; ^ 00 , and lead to instabilities. Except when otherwise stated, the analysis in this section will hold for 
general values of to^, to^ and A. 
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A. The near region solution 

In the near-region, k^x <C |1 — one can neglect the k^x term, and the radial equation (16) is approximated by 


x{l + x)d‘lh + (1 + 2x)dxh + - 


1 - 


^2 ^2 


a; -I- 1 X 


h = 0. 


(71) 


With the definition h = a;l‘”l/^(l-|-a;)^/^ w, the near-region radial equation becomes a standard hypergeometric equation 
[39] of the form 


a;(l -I- x)d'^w + [c+ (a -I- 6 -I- l)x]dxW + abw = 0, 


(72) 


where 


a — - (1-I-Id + ^ + J^) , ^ — 2 + 1^1 + ^ ~ > c—1-1-Id- (73) 

The full solution to the above is given in terms of hypergeometric functions as w = AF{a, b, c, —x) + B x^~‘^F{a — 
c + l,b — c + 1,2 — c, —x), which allows us finally to write the solution of the radial equation in the near region as 

h = Ax^‘’^^‘^{l + x)^^‘^F{a,b,c,—x) 

+B a;“l‘’l/^(l -I- x)^^‘^F{a — c-|-l,6 — c-|-l,2 — c, —x). (74) 


At this point we impose the first boundary condition: the solution must be regular at a; = 0 since the geometry is 
smooth at the origin of the “core”. The term proportional to diverges at a; = 0 and must be discarded, i.e., 

its coefficient, B, must be set to zero. 

To perform the matching we need to know the large x behavior behavior of the regular near-region solution. To 
this end, one uses the x —>■ 1/x transformation law for the hypergeometric function [39] 

F{a,b,c,-x) = —^a:"“F(a, 1-c-ka, l- 6 -ka,-l/a;) 

i (oji (c — a) 

+ — ^x~^ F{b,l-c+b,l-a+b,-l/x), (75) 

1 (a)l (c — o) 


and the property F{a,b, c,0) = 1. Note that this expression for the transformation is only valid when a — b = v is 
non-integer. This is an assumption we will continue to make throughout this section. In the end, we shall derive a 
condition determining the allowed eigenvalues that will not be dependent upon this assumption and therefore we may 
extend our results to integer values of v by continuity. 

The large x behavior of the near-region solution is then given by 


h 


Ar(i + |d) 


■_THd_ 

.r[Hi + ICI + ?-^^)]r[i(i + |d-e-^^)] 

,__ 

r[Hi + ICI + C + ^)]r[Hi + ICI-e + ^)] 


(76) 


B. The far region solution 


In the far-region, kx"^ ^ max{^^ 
can be approximated by 


1 , the terms ^^/(x -|- 1 ) and (“^/x can be neglected, and the radial equation 


d^ixh) 



ly^-l 

4x2 


(xh) = 0 . 


(77) 


The most general solution of this equation when v is non-integer is a linear combination of Bessel functions of the 
first kind [39], 


h = x [CJi,{k\/x) + DJ-,^{k\/x)] . 


(78) 
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This form does not lend itself easily to application of the boundary conditions. Instead, for large k^/x, the solution 
may be expanded as [39] 


CC-3/4 

^/2TrK 


-i 


e (^Ce *”2 + Ilex’s ) + 




De- 


(79) 


As in the WKB analysis, we assume that the real part of oj is negative, and therefore the positive and negative sign 
exponentials give, respectively, ingoing and outgoing waves. We require that there be purely outgoing waves at infinity 
and so impose the constraint that the coefficient of the positive exponential vanishes, yielding 


C = 


(80) 


When uj becomes complex, so too does k. Since the sign of the real part of oj is negative, the dehnition of k (17) 
implies that its imaginary part has a sign opposite that of the imaginary part of uj. Therefore, requiring additionally 
that the solution be finite as a; ^ 00 implies that the imaginary part of u must be positive. This is precisely the sign 
for the imaginary part of the frequency that leads to instabilities. Thus we see that simply requiring the solutions 
with complex frequency be finite at infinity automatically guarantees they lead to instabilities. 

Now, to do the matching in the overlapping region, we will need to know how the far-region solution behaves for 
small values of x. More specifically, for small k^/x, and considering only the dominant terms, the solution behaves as 
[39] 


h 


iv^r 

r(i + j.) 


1 / 


u-l 

X 2 


r(i - u) 



(81) 


C. Matching conditions: Selection of frequencies 


We will now determine the frequencies that can appear when the geometry is perturbed by a scalar field. The 
frequency spectrum is not arbitrary: only those values that satisfy the matching conditions between the near-region 
and the far-region are allowed. We shall see that there are two solutions of the matching equations, yet only one will 
lead to instabilities. 

Matching the powers of x between the near (76) and far-region solutions (81), and taking a ratio to eliminate the 
amplitudes A and D, yields 


2. r(i-i2) ^ r(t2) r(i(i - ^ + |CI + 0 ) r(i(i -+ |ci - g)) 
r(i + r(-^) r(i(i + ^ + |CI + ?)) r(i(i ++ |c| - 0 ) ' 


The problem of hnding the outgoing modes thus boils down to solving the single transcendental equation (82); we will 
do so by iteration. Note that the k, dependence on the left hand side means that it is suppressed. For the equation to 
hold, a similar suppression must also occur on the right hand side. This is only possible if one of the gamma functions 
in the denominator of the right side is large. Since the gamma function diverges when its argument is a non-positive 
integer, we take as a first iteration the choice 


^+IC|-C = -(2iV + l) , (83) 

where the non-negative integer N will again be referred to as the harmonic. Note that we could also have chosen 
the above relation, but with the opposite sign for While this does indeed lead to a solution, one finds that the 
imaginary part of the frequency is always negative, i.e., the modes are exponentially damped in time. 

This first estimate is obviously not the end of the story as it would cause the right side to completely vanish. To 
go beyond this approximation, we rewrite Eq. (82) in terms of N, then perturb N ^ N + SN, where SN <C N. 
This deformation appears at leading order only for the F function in the denominator on the right hand side that 
diverges, it may be neglected in all other factors. More concretely, to extract SN from the F function we use 
F(z)r(l — z) = 7r/sin(7r2), and sine function identities to obtain the expansion 

r{-N - SN) -[{-1)^N\SN]~\ (84) 

Substituting this into (82), and using a number of F function identities, we solve for the imaginary part of the first 
correction 

lm{SN) = V'U Hw-eici > 


(85) 
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where Mn=n”=l(l + a/i). Since N is 0(1) and 5N ~ we see that we may stop after the first iteration. As a 
function of v, this can have a single maximum near v ^ k. In general we will have k <C 1 and ~ 1 + /, so will always 
be in a region where this is a monotonically decreasing function of v. For fixed v, the last two factors make this an 
increasing function of N and ICj, but the general behavior will be dominated by the effects of changing v. 

The equation (83) uniquely determining ui can be exactly solved 


+ Q^)uj 


f , Cp ^ 



f \ ^p^p^ 

n'— 



2 






( 86 ) 


where 




+ M(Si + Sg + Cp)) 


c = eo-|C|-(2iV + l) , 


(87) 


and a variable with a subscripted 0 means we have set oj = 0. Note that as long as m > n + 2, we have s/g"^ <C 1 
and both quantities are positive. When m ^ n + 1, though, e —oo (since M ^ 0, ^ —oo and stays 

finite), ensuring that there can be no instability for the supersymmetric solutions. This extends to arbitrary modes 
the conclusion from the discussion associated to equation (40) for modes with = A = 0. 

When evaluated on a solution, v is given hy v = ug + c. Since we are interested in solutions for which ui is negative, 
this means c > 0. Then, requiring that lo be negative and real, gives three more conditions. The first ensures that the 
result is real while the second requires that the first term of (86) is negative. Finally, the condition that appears to 
be the most difficult to satisfy ensures the contribution from the square root does not make the total result positive, 
i.e., 


-i'n>0 


( 88 ) 


When A ^ 0, these conditions must also be supplemented by the requirement that — A^ > 0, which ensures the 
asymptotic behavior of the solution is correct. With these satisfied, we may determine the effect of the correction. 

The imaginary contribution to N is taken as resulting from a small imaginary correction to u. Then, the two are 
related through 


6 N 


Sivd 
Scv 


2iy 


N 


H“ s)uj -\- {XspCpM/-\- ^c)] 




(89) 


In the final line we have used the solution (86) to show that the sign of 5N determines the sign of the correction to 
to. Since Jm{6N) is always positive when evaluated on the solution of (83), the corresponding imaginary part of to is 
positive. 

To summarize, whenever the constraints, in particular (88), are satisfied there is a corresponding outgoing mode 
of the scalar field equation. Further, the imaginary part of the frequency of this mode is guaranteed to be positive, 
indicating that it leads to an instability. The timescale for the instability generated by the mode is a monotonically 
increasing function of v, which is given by 


/ 2n X M 

{e+g)i> = ec- XgSpCp^ + 


^(^ec - XgSpCp^^ + (£ + g'^){2cgXspCp^ + v^g^ - ec^) . 


(90) 


A similar argument, based on the solution of equation (83), but with the opposite sign for ^ would lead to a set of 
outgoing modes with an amplitude that decays in time. 

As an example, consider the particular background geometry and scalar field solution described by 


m = 5; n = l; ai = 19.1; ci = 5 ; Cp = 1.05 ; 

A = rUfj, = 0; I = = 2 . (91) 

The first two iterations with N = 0 gives 


uj = -2.8717, 

T~^ =Tm{5uj) = 4.42 x 10“^^ , 


(92) 
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The results obtained here are consistent with the WKB analysis of the last section, i.e., there are outgoing modes 
that rotate in the same sense as the background geometry whose amplitude grows exponentially in time. What we 
have gained is an explicit set of relations that allows the unstable mode frequencies to be calculated. In particular, 
one can now make definite statements about the relative timescales for unstable modes just by looking at equation 
(90). We leave the precise details of this to Appendix C and just give the results here. The most unstable modes are 
those which minimize v. Since e this generally means that the modes which maximize c or minimize vq will be 
the most unstable. In general this means we should consider the lowest possible I for which the constraints can be 
satisfied when setting my, = I, my, = 0 and N = 0. 

A second benefit of this analysis is an improvement in accuracy for the most unstable modes. For comparison, 
performing the WKB analysis and not neglecting any terms in the potentials or approximating the bottom of the well 
with a parabola gives oj = —3.129 + 4.00 x 10“^°z. From the full numerical solution we have to = 2.8718 + 4.46 x 10“^^z. 
For values of uj in this range we have k~'^ ~ 1900, so we are well within the range for which we should trust this 
solution. As k~^ approaches max(|^|,,J), this analysis begins to break down, but it appears that the WKB approach 
becomes increasingly accurate. In the next section we will present a more detailed list of eigenvalues corresponding 
to instabilities and discuss the results. 


VI. NUMERICAL RESULTS 

We will now solve the radial equation (16) numerically to extract the instability. We begin with an exposition of the 
numerical algorithm. The only approximation used in this section concerns the angular eigenvalue A, that we shall 
assume to be well described by (19). At the end of the calculation we always make sure the result fits the regime of 
validity of this approximation. Note, however, solutions can still be found even when outside this range. The easiest 
way to do this is by treating A and to respectively as eigenvalues of the angular and radial equations. The coupled 
system may then be solved by first assuming the approximation to hold and solving the radial equation for uj, this is 
then fed into the angular equation to obtain an improved value of A. This process may be iterated until the desired 
level of convergence is achieved. 


A. Numerical procedure 

The method of finding solutions numerically is very much like performing the matched expansions. We use Eqs. (74) 
and (78) to Hx the initial conditions for two integrations of the exact radial equation. Since the equation of motion 
is linear, we may immediately match the two solutions at a point in the interior region by rescaling. This leaves two 
more conditions to be satisfied, that matching of the derivatives of the real and imaginary parts. Fixing all other 
parameters, we vary the real and imaginary parts of oj to satisfy these conditions. 

Given the small size of the expected imaginary part, it is most straightforward to use a package like Mathematica 
[40] with its software based arbitrary precision, to perform the calculations. Satisfying the matching conditions can 
be done by treating the difference in derivatives at the interior point as a complex valued function of w. A root may 
then be searched for using the built-in function FindRoot which, for a function without explicit derivatives, looks for 
the solution by constructing secants. 

Since the imaginary part is expected to be far smaller than the real, gradients of the matching function in the 
imaginary oj direction will be large only when very near a solution, but negligible elsewhere. The initial guesses at 
the solution are therefore very important for ensuring that iterations converge to a solution. It was found empirically 
that solutions could consistently be found by choosing to start the search in a region around the real value of oj for 
which the inner solution vanishes at the matching point. Small changes in the imaginary part of oj near this point 
appear to be sufficient to bring about convergence. In Fig. 2 we show an example solution obtained in this manner. 
The solid line is the full numeric solution, with the integration starting at small x in red to the left of the black dot 
and that starting at large x on the right in blue. The dashed lines are the near (74) and far (78) approximations 
used to set the initial conditions for integrating the exact radial equation. The fact that the imaginary part of oj is in 
general very small raises non-trivial problems, related to the number of digits of precision used and the exact way in 
which boundary conditions are applied. A discussion of these aspects is deferred to Appendix B. 
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FIG. 2: An example solution showing vanishing as both a; ^ 0 and x oo. 


B. Numerical results 

Our numerical results are summarized in Fig. 3 and Table I. In Fig. 3 on the left we present the numerical solutions 
obtained for 


m = 5 , n = 1, Cl = 1.1, C 5 = 1.52 , oi = 262.7, A = = 0 . (93) 

where we consider only the lowest harmonic, N = 0, but vary I = m^. At ^ = 1, ^ 40 indicating the matched 

solution is valid, as I a grows so do ^ while k~^ shrinks, meaning the approximation should soon break down. At 
I = 5, K~^ ~ 10 and the approximation is becoming no longer valid. Finally, when I = 13, k~^ ~ 1 and differences 
between the matched and numerically determined eigenvalues are starting to become apparent. In Fig. 3 on the right, 
we use the same parameters as before, but now fix I = = 4 and vary the harmonic from iV = 0 up to 4. Increasing 

N leads to smaller values of w and therefore smaller values of k, so that the matched solutions are valid throughout. 
It should also be noted that if the approximation a\uP/B? <C is valid for a given then it should be valid for 
all m^. This is because uj scales with m,/,, as we observed within the WKB approximation. 
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FIG. 3: On the left we choose the lowest harmonic and vary I = from 2 to 13 from upper right to lower left. The solid circles 
represent the numeric solutions, while the triangles are the results of the WKB analysis and the unfilled circles correspond to 
the matched expansion. On the right we fix 1 = 4 and vary the harmonic from 0 to 4 from upper left to lower right. 

In Table I we present and compare the numerical results with those obtained through the approximate analyical 
approaches. The values labeled as WKBnum (numerical WKB) stand for values obtained using the full WKB 
approximation, formulae (48), (49), (51) and (55), which have been handled numerically. The values obtained using 
the parabolic approximation, formulae (56)-(59), for the potential are denoted by WKBan- 
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Numeric 

WKBnum 

WKBan 

Matching 

1 

-0.184-b 3.83 X 10“®i 

- 

- 

-0.184 + 3.83 X 10"®i 

2 

-0.744-b 2.51 X 10“®i 

-0.812 + 1.61 X 10“’’i 

-0.826 + 1.89 X 10~'^i 

-0.744 + 2.64 X 10"®i 

3 

-1.312-b 3.73 X 10-^i 

-1.359 + 1.20 X 10“®i 

-1.371 + 1.48 X 10"®i 

-1.312 + 3.53 X 10"®i 

4 

-1.883 -b 3.69 X 10~i°i 

-1.919 + 9.33 X 10-“i 

-1.932 + 1.17 X 10"®i 

-1.882 + 3.63 X 10-i°i 

5 

-2.456 -b 3.55 x 10"“i 

-2.486 + 7.37 x 10”“i 

-2.499 + 9.39 x 10'“i 

-2.454 + 3.39 x lO’^i 

6 

-3.030 + 3.22 X 10-^^i 

-3.055 + 5.89 X lO’^^i 

-3.072 + 7.62 X lO'^^i 

-3.028 + 3.02 X 10-^^i 

7 

-3.605 + 2.77 X lO’^’^i 

-3.626 + 4.73 x lO’^’^i 

-3.647 + 6.23 x lO’^’^i 

-3.602 + 2.63 X 10"^®i 

8 

-4.180 + 2.47 X lO’^'^i 

-4.199 + 3.82 X lO’^'^i 

-4.216 + 4.88 X lO-^'^i 

-4.176 + 2.24 X 10~^'‘i 

9 

-4.755 + 2.05 X 10~^®i 

-4.772 + 3.09 X 10"^®i 

-4.794 + 4.03 X 10"^®i 

-4.751 + 1.89 X 10"^®i 

10 

-5.331 + 1.76 X 

-5.346 + 2.51 X 10”^®i 

-5.369 + 3.26 x 10'^®i 

-5.326 + 1.58 X 10"^®i 

11 

-5.907 + 1.49 X lO’^'^i 

-5.921 + 2.03 X lO’^'^i 

-5.947 + 2.65 x 

-5.902 + 1.32 X 10"^'^i 

12 

-6.483 + 1.22 X 10"“i 

-6.496 + 1.65 X 10"“i 

-6.516 + 2.07 X 10"^®i 

-6.477 + 1.09 X 10-^®i 

13 

-7.059 + 1.04 X 10-^®i 

-7.071 + 1.34 X 

-7.102 + 1.81 X 10"^®i 

-7.053 + 8.97 X 10"^“i 


TABLE I: Some numerical values of the instability for a geometry with m = 5, n = 1, ci = 1.10, C5 = 1.52, ai = 262.7, A = 
m^, = 0 and I = m,/,. In the second column, we have the results of the full numerical analysis; in the third column, WKBnum 
(numerical WKB) stands for values obtained using the full WKB approximation, formulae (48), (49), (51) and (55); and in the 
fourth column, labelled as WKBan, the values obtained using the parabolic approximation for the potential, formulae (56)-(59), 
are given. In the final column, we present the results of the matching procedure (83),(85). Notice the close agreement between 
all the different methods. For m,/, = 1 — 1 and for these particular values of the parameters, the WKB analysis, as done here, 
breaks down. Indeed, for = 1, the potential V+ has no minimum. 


Notice first that all the different approaches yield consistent and in fact very similar results: they are all rather 
accurate in their own regime of validity. As predicted by the analytic approaches, and verified numerically, the real 
part of the frequency scales with m^, whereas the logarithm of the imaginary part scales with e.g., see Eq. (49). 
Thus the instability timescale increases rapidly as a function of 


VII. DISCUSSION 

In this paper, we have shown that the non-supersymmetric JMaRT solitons [1] are classically unstable. The relevant 
instabilities are quite generic to spacetimes which have an ergoregion but are horizon-free [26]. However, as noted in 
Section IIB, the general proof does not strictly apply to the JMaRT solutions since the latter support nonradiative 
negative energy modes as shown in Appendix D. Hence we have explicitly shown that the ergoregion instabilities are 
active in the JMaRT geometries using three different approaches, which in the end show a remarkable agreement — 
see Fig. 3 and Table I. Perhaps the most physically intuitive method is the WKB analysis carried on in Sec. IV. This 
approach allows us to clearly identify the nature and physical properties of the instability. However, this analysis 
is only expected to be valid for large angular momentum quantum numbers, i.e., rriTp ^ 0, which is not where the 
instability is strongest. The more unstable modes were studied using the matched asymptotic expansion method [37] 
in Sec. V. As a final consistency check of these analytical results, we made a numerical analysis of the wave equation 
in Sec. VI. 

In passing we note by considering orbifolds, the JMaRT solutions were extended to a six-parameter family which 
includes a third integer k characterizing the orbifold group [1]. Of course, it is straightforward to adapt our 
instability analysis so that the modes respect this orbifold symmetry in the covering space and so one concludes that 
the ergoregion instability arises in these orbifold geometries as well. 

Let us now summarize some of the features of the ergoregion instability found for the JMaRT solutions: 

(i) The general shape of the WKB potentials V± are sketched in Fig. 1 for the case in which an instability is present. 
The key point is that when the ergoregion is present the bottom of the potential well in V+ reaches negative values. 
The unstable modes are those whose pattern speed is negative and approaches the minimum of V+ from above 
(see Fig. 1). Thus, they are nearly bound states of the potential well in V+ that can however tunnel out to infinity 
through V-. 

(ii) The fact that the unstable modes are those with negative phase velocity, < 0, has a clear physical interpretation. 
As in the discussion of Eqs. (32) and (33), modes with E^ < 0 are those that propagate in the same sense as geometry’s 
rotation Therefore at infinity these modes carry positive angular momentum (same sense as 12^), as well as positive 
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FIG. 4: Damped modes are those that have positive 


energy. Hence by conservation of energy and angular momentum, with the onset of the ergoregion instability, the 
JMaRT solutions are shedding both energy and angular momentum by an amount that increases exponentially. 

(iii) The instability can be quite strong, depending on the particular combination of parameters that define the 
geometry. More importantly, the instability is robust, in the sense that it exists for a wide range of parameters. 

(iv) With m = n + 1, the JMaRT solutions are supersymmetric and so must be stable. It is a consistency check of our 
analysis then that we hnd no instability in this case. As commented in Section IV, when m = n + 1 the potential V+, 
as given by Eq. (40), is always positive. Hence there are no negative modes which could intersect the potential 
well of V+ and the SUSY geometry is stable as required. 

In our analysis, we have focused on the special case A = 0 and = 0, to simplify the relevant equations. In 
fact, the ergoregion instability persists when either or both of these parameters are nonvanishing. A discussion of the 
general situation is given in Appendix C. The result is most simply understood from the point of view of the WKB 
approach. Then all of the additional contributions to the effective potential (29) introduced by a nonvanishing or 
A are suppressed by inverse powers of and so can certainly be neglected in the limit of large One can further 
check that the instability exists over some range even when does not dominate the other two. One distinguishing 
feature of A 7 ^ 0 is that asymptotically the scalar modes have an effective mass in five dimensions. In our analysis, 
this is reflected in the fact that asymptotically V± —> ±|A/m^| and so there is an additional barrier for the modes 
to tunnel out to infinity. However, for sufficiently large m^, such tunnelling is possible. One other interesting point 
about the large regime is that unstable modes appear with either sign of and A. Hence while for the modes 
on which we have focussed, the instability is ‘powered’ by and results in decreasing this angular momentum, there 
are unstable modes which may at the same time increase \ J^\ and/or P. 

Adding a mass for the scalar held modihes the potentials V± in essentially the same way as having nonvanishing 
A. Hence we expect the ergoregion instability will even appear for massive helds, at least in modes with sufhciently 
large angular momentum. As described in section HB, the arguments given by Friedmann [26] are quite general and 
so we expect the ergoregion instability to appear for higher spin helds as well. In particular, we expect the helds 
of the low energy type Hb supergravity will generically experience this instability. Having said that the ergoregion 
instability is robust, we must also add that it can be suppressed in certain parameter regions. In particular, one hnds 
that the instability timescale becomes extremely long in the regime where Qi and Q 5 are much larger than the other 
scales. Further we add that in the decoupling limit where one isolates an asymptotically AdSa core [1], the ergoregion 
instability is absent. The simplest way to understand this result is that the AdSa core has a globally timelike Killing 
vector [1] and so there is a ‘rotating’ frame where we can dehne all energies to be positive. One can also explicitly 
verify the absence of an ergoregion instability in the core solutions by directly applying the analysis used in this paper 
to those backgrounds. 

The JMaRT geometries (both supersymmetric and non-supersymmetric) also have damped modes, i.e., modes (12) 
for which the imaginary part of oj is negative. As per the WKB analysis, these are modes with positive below the 
local maximum of U+ that tunnel out to infinity through V+ — see Fig. 4. 

As emphasized previously, we can also find purely bound states {i.e., nonradiative modes) with oc — A^ < 0. 
With some fine-tuning, it may also be possible to find geometries which support bound states with n = 0. These 
nonradiative modes are described in Appendix D. The typical situation for such modes is sketched in Fig. 5. As 
already noted above when A 7 ^ 0, asymptotically V± ±|A/m^| and so there is a finite potential barrier at infinity. 
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FIG. 5: Qualitative shape of the potentials V+ and V- when < 0. These are the purely bound states that are discussed 

in Appendix D. 


If this barrier is sufficiently large relative to = w/m^, bound states can arise. These bound states can also be 
negative energy states, as can be seen with the energy integral (7). The absence of such negative energy modes which 
do not radiate at infinity was central to Friedman’s general argument for the ergoregion instability. In [26], he did 
not find any such nonradiative modes because he only considered the massless fields for which there is no potential 
barrier at infinity. Note, however that the current situation is more complicated because the KK-momentum of the 
background, as well as the angular momenta, contribute to the presence of an ergoregion. 

The appearance of negative energy states in the presence of an ergoregion can be anticipated from a geodesic 
analysis [27]. By definition, the Killing vector which generates asymptotic time translations, becomes space-like 
inside the ergosphere. Hence (time-like or null) geodesics can have either positive or negative energy, e = —t ■ u, in 
this region. However, asymptotically only positive energy (i.e., future-oriented) geodesics are physical. Therefore any 
negative energy geodesics must be confined to circulate within the ergoregion. Of course, in a black hole background, 
such geodesics would ‘disappear’ behind the event horizon. However, for horizon-free geometries, such as the JMaRT 
solutions, they are stable bound orbits and so it is natural to find bound states in the context of a field theory analysis. 
However, the question then becomes whether the analogous modes of the field ‘fit’ inside the ergoregion or whether 
they ‘leak’ out to infinity, i.e., whether a negative energy bound state or an ergoregion instability results. A more 
thorough examination of the bound states shows that the negative energy bound states states are characterized by 
having Ej, = w/A < 0 while the ergoregion instability is associated with modes where E^ = and/or E^ = 

are negative - see Appendices C and D. Hence as the geodesic analysis would suggest the negative energy modes 
have a negative pattern speed or phase velocity, but the KK-momentum modes tend to lead to bound states while 
the spinning modes are related to instabilities. 

The presence of negative energy bound states can also be expected to enhance the decay of these horizon-free 
geometries. The analysis of the ergoregion instability (considered in this paper) is only at the level of linearized 
test fields. Generically any theory coupling to gravity will also have nonlinear interactions (e.g., even the free scalar 
considered here has nonlinear couplings with gravitons). These nonlinear couplings might be expected to lead to 
processes, where positive energy modes are radiated at infinity while negative energy modes are populated within 
the ergoregion. However, one should note that the negative energy modes are exponentially decaying at large radius 
— see Appendix D — while the positive energy modes are power-law suppressed inside the ergoregion. Hence the 
overlap of these modes is expected to be small, which will suppress this nonlinear contribution to the decay. 

We now turn to consider the endpoint of the ergoregion instabilities. As emphasized before, the presence of these 
instabilities relies on two key ingredients, namely, the geometry has an ergoregion but it does not have an event 
horizon. Hence the resulting decay process could be terminated either by the disappearance of the ergoregion or the 
appearance of a horizon. However, the unstable modes radiate with a positive energy density asymptotically which 
is compensated for by a negative energy density inside the ergoregion — as could be seen in Eq. (7). Hence the 
onset of the ergoregion instability produces a(n exponential) build-up of negative energy near the core of the JMaRT 
solutions. Therefore it seems unlikely that an event horizon will form since the latter is typically associated with 
a large build-up of positive energy density. This reasoning then suggests that the decay must terminate with the 
disappearance of the ergoregion. The supersymmetric D1-D5-P microstate geometries [15, 18-20] are all free of an 
ergoregion and hence it is natural to suppose that these are at the endpoint of the ergoregion instabilities. Of course. 
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these solutions offer a huge family of possible endpoints and the precise one that forms will depend on the details of the 
decay process, beyond the linear regime considered here — although as we are only considering the classical evolution, 
it is in principle possible given a certain set of initial conditions. Of course, we can expect that the final mass should 
be close to the BPS mass determined by the charges of the initial JMaRT solution, i.e., E = Tr/dGs [Qi + -\- Qp\. 

Although even here, we can only say ‘close’ as we know that the unstable modes with A 0 (and either sign of 
A) occur which may modify the final value of Qp. Similar comments apply for the angular momenta, J^p and Jp,. 
We also observe that there is no reason to expect that the decay process will lead to an endpoint within the family 
of supersymmetric JMaRT solutions. Of course, at the level of the present discussion, we cannot rule out that the 
endpoint is only a nearly supersymmetric solution (or that this would be the effective endpoint). Our expectation is 
that such solutions will have a ‘small’ ergoregion and that the instability might be eliminated (or strongly suppressed) 
before the ergoregion precisely vanishes. 

The stability analysis of the JMaRT solitons [1] is relevant for the stringy tachyon decays discussed recently in [29]. 
Originally, [30] considered tachyon condensation in certain D1-D5 black string backgrounds where tachyonic string 
winding modes can occur if one chooses antiperiodic boundary conditions for the fermions around the circle on which 
the black string is compactified. The latter choice necessarily restricts the scenario to a non-supersymmetric sector 
of string theory which already suffers from various instabilities [41]. Ref. [29] considered adding angular momentum 
to the black strings. In this case, it was shown that string modes winding certain compact circles near the horizon 
can be tachyonic even when the asymptotic fermion boundary conditions are supersymmetric. The relevant point for 
the present discussion is that the endpoint of the tachyon condensation is in general one of the non-supersymmetric 
JMaRT solitons. Now, in this paper, we have shown that these solitons are themselves unstable and so they will 
not be the final endpoint of these decays. Instead, the ergoregion instability will continue the decay process and as 
suggested above, will likely terminate with a supersymmetric microstate geometry. 

We would now like to consider the implications of ergoregion instabilities for Mathur’s fuzzball program of describing 
black holes in terms of an ensemble of microstate geometries. If this program is to succeed it must supply a description 
of both supersymmetric and also non-supersymmetric black holes. At first sight, it may appear that constructing 
non-BPS microstate geometries is not possible. In particular, non-BPS states will decay and so it is not clear that 
there should be stationary geometries to describe them. However, the JMaRT solutions provide an explicit example 
indicating that this is not really a problem. In fact, the decay of non-BPS microstates was already considered in the 
D-brane description of nonextremal black holes [3]. In that context, it was seen as a success of the string theoretic 
approach as this instability had an interpretation in terms of Hawking radiation [42-44]. Of course. Hawking radiation 
is a quantum effect in the black hole background and so presents no obstacle to the construction of classical supergravity 
solutions which are static or stationary. 

It is perhaps useful to remind ourselves as to how this distinction arises. The classical limit can be understood 
as the limit where the string coupling is vanishingly small [43]. However, the interesting classical solutions are 
those which correspond to states where the various quantum numbers are extremely large. That is, ni, ns oc I /gs 
and rip, oc f/gl while gs —> 0. These scalings are chosen to ensure that the gravitational ‘footprint’, i.e., Qi, 

Qpi 0,1 and 02 , associated with each of these quantum numbers remains finite. However, in this limit, the ADM 
energy of the system diverges with E oc f/g^- As the energy is a dimensionful quantity, this can be accommodated by 
changing the units with which energies are measured in the classical limit. Essentially this divergence is associated 
with the divergence of the Planck mass, which does not serve as a useful reference scale in classical gravity. Now the 
decay rate of the nonextremal DI-D5-P black holes can be computed in a straightforward manner [43, 44]. The key 
point, however, is that the final expression for dE/dt is expressed in terms of geometric quantities and is independent 
of gs- Therefore in the classical limit, the rate of energy loss remains finite in stringy units but becomes vanishingly 
small when measured against the fiducial energy scale that was established for classical physics. Hence the non-BPS 
black holes become stable in the classical regime. 

We note that the ‘straightforward’ calculations of the decay rate referred to above can be performed either in 
the framework of a microscopic D-brane perspective or of the gravitational perspective of Hawking radiation. The 
suprising result is that the results of both analyses agree precisely [43, 44], including greybody factors, at least in the 
so-called ‘dilute gas’ approximation [45]. However, even though suggestive arguments can be made in this regime [46], 
this remarkable agreement remains poorly understood. As the JMaRT solutions are horizon-free, the gravitational 
calculation of the decay rate would have to be modified. Using the connection between absorption and emission 
rates, it is possible that absorption calculations along the lines of those presented in [I] could be extended to yield 


We should note that the JMaRT solutions begin in a low-mass regime where < (ai — 02 )^, however, if the ergoregion instability 
sheds the background angular momentum efficiently then the system will evolve to a regime where black holes can form. Hence we can 
not rule out the appearance of an event horizon - we thank Simon Ross for correspondence on this point. 
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the desired decay rate. On the other hand, the underlying microscopic states for the JMaRT solutions were already 
identified in [1]. Hence one can use microscopic techniques to estimate the decay rate expected for these solutions. 
The result is dE/dt ~ QiQ 5 {m — n)^/R^ and again this quantity remains finite as gs 0. Therefore we can again 
ignore this decay channel for the classical JMaRT solutions. 

However, the ergoregion instability investigated in this paper is a classical instability and so should not be associated 
with the decay discussed above. We should also note that the form of these two instabilities differs. Above one is 
considering the spontaneous decay of the system while the classical instability really corresponds to a decay that 
results when the initial data does not precisely match that of the JMaRT solutions. Of course, in the quantum 
regime, the same modes associated with the ergoregion instability will give rise to spontaneous decay due to quantum 
fluctuations of the background.However, the latter will again be suppressed in the gs ^ 0 limit. This reflects the 
fact that the background can be prepared with arbitrarily accurate precision in the classical limit and so it should be 
possible to produce an arbitrary suppression of ergoregion instability. Alternatively, working in the classical limit, we 
can regard the ergoregion instability as a property of how the JMaRT solutions interact with external sources. That 
is, generically if an external wave packet impinges on one of the non-supersymmetric JMaRT configurations, it will 
produce a dramatic decay of the original background. Hence this instability seems to present a major challenge for 
the fuzzball description of black holes. 

We have argued that the ergoregion instability is a robust feature of the non-supersymmetric JMaRT solutions 
over a wide range of parameters. Given general arguments along the lines of [26], we also expect that this instability 
will be a generic feature of any smooth horizon-free geometries which describe microstates which are non-BPS and 
carry significant angular momentum (and hence have a macroscopic ergoregion). Therefore if a nonextremal D1-D5-P 
black hole is to be described by a coarse-grained ensemble fuzzball, it seems that that the classical black holes must 
suffer from an analogous instability. While the presence of an event horizon eliminates the possibility of an explicit 
ergoregion instability, there are, in fact, a number of potential instabilities which might afflict these black holes and 
possibly reproduce the same physics: 

a) Superradiant Instability: Spinning nonextremal black holes will exhibit superradiant scattering, where an 
incident wave packet can be reflected with a stronger amplitude. Superradiance by itself does not provide a classical 
instability, but an instability can arise if the scattered modes are reflected back to rescatter, as described in Section 
H. This scattering was considered for higher dimensional spinning black branes [33] and there it was found that when 
the noncompact space has more than four dimensions, this instability does not arise. Explicitly analyzing the present 
D1-D5-P black string again seems to indicate the absence of an instability [47]. 

b) Gyration Instability: Considering supersymmetric D1-D5-P black strings, it was found that above a certain 
critical angular momentum a straight black string is unstable towards carrying the angular momentum in gyrations 
of the horizon [48]. This instability should also appear in non-supersymmetric configurations and so would present 
an instability at large values of the angular momentum. 

c) Gregory-Laflamme Instability: The relevant configurations are black strings and so are expected to suffer from 
the Gregory-Laflamme instability [49] in two ways. The first is the usual instability of long wavelength modes along 
the string. Of course, this instability can be eliminated by reducing the radius of the compactification along the string. 
For a fixed radius, it is also suppressed by the boosting along this direction which induces the KK-momentum [50]. 
This instability is not related to the angular momentum carried by the black string or the presence of an ergoregion, 
but we list it here for completeness. 

d) Ultra-spin Instability: In six or higher spacetime dimensions, one can find black hole solutions with an arbitrarily 
large spin per unit mass [36]. However, it was argued [51] that a Gregory-Laflamme-like instability will arise to 
dynamically enforce a Kerr-like bound in these cases. While this analysis does not directly apply in five dimensions, 
entropy arguments suggest an analogous instability still exists and will lead to the formation of a black ring if the 
angular momentum is too large [52]. 

While there are several possibilities for instabilities of a black string in six dimensions, it seems that none of these 
can reproduce the physics of the ergoregion instability which will afflict the non-BPS microstate geometries. This 
observation relies on the fact that these instabilities have a different character at a very basic level. The ergoregion 
instability might be termed a radiative instability, in that, the instability is by definition connected to modes that 
radiate at infinity. In contrast, the four instabilities considered above for black strings can be termed internal 
instabilities. That is, these instabilities are primarily associated with a rearrangement of the internal or near-horizon 
structure of the black string. While these instabilities will be accompanied with some radiation at infinity, this 
will be a secondary effect with these instabilities. Therefore it seems that emulating the ergoregion instability in a 
nonextremal black string background will require the discovery of a new kind of instability. While we are performing 


In [1] it was erroneously assumed that all of these geometries have an AdSs core to argue that such emissions would not occur. 
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a detailed analysis of the nonextremal D1-D5-P black string, our preliminary results indicate that no such instability 
arises [47]. 

We also note in passing that at the same time the microstate geometries should be able to emulate any instabilities 
found in the black string backgrounds. In particular, the Gregory-Laflamme instability is a robust instability that will 
afflict these backgrounds for sufficiently large R. In the microstate geometries, one should then find unstable modes 
carrying KK-momentum which are confined near the core of the soliton. We have studied bound states for a test 
field in the JMaRT solutions, as described in Appendix D. While the modes we identified only arise for nonvanishing 
KK-momentum as desired, they are all stable, i.e., they have real frequencies. Hence they can not serve as the analog 
of the Gregory-Laflamme instability in the non-supersymmetric JMaRT solutions. However, the latter would be a 
gravitational instability, i.e., it should not be expected to appear as a scalar test field, and so this question requires 
further investigation. 

A possible reconciliation of these ideas with the fuzzball proposal would be that the microstate geometries could 
provide an accurate description of a black hole but only over a long but finite time. In the context of the AdS 3 /GFT 2 
duality, some evidence for such a picture has recently been found [53]. With this new point of view, a key question is 
to determine the timescale over which microstate geometries cannot be distinguished from black holes. One suggestion 
[53] is that it should be of the order of the recurrence time, which would be exponential in the relevant quantum 
numbers. An alternative suggestion might be that the timescale is associated with Hawking evaporation which would 
involve (inverse) powers of the quantum numbers. However, note that both of these suggestions diverge in the classical 
limit. Hence the ergoregion instability found here seems to be in conflict with both of these suggestions. While the 
instability timescale is certainly very long in certain parameter regimes, it is a classical timescale, i.e., it is finite in 
the classical limit. Hence our results would suggest that spinning microstate geometries and black holes should be 
distinguishable on a large but classically finite timescale. 

However, one must ask how characteristic our results for the JMaRT solutions will be of generic microstate ge¬ 
ometries. In particular, we note that the GFT states corresponding to the JMaRT solutions are exclusively in the 
untwisted sector [1, 19, 54]. On the other hand, the majority of microstates accounting for the entropy of the black 
strings are expected to be in a (maximally) twisted sector [11]. From a geometric point of view, we would observe that 
the JMaRT solutions have all the same Killing symmetries as the D1-D5-P black holes, while the generic microstate 
geometry is expected to have a complex nonsymmetric core. Therefore it is not unreasonable to expect that the 
ergoregion instability timescales found for the JMaRT solutions will not be characteristic of the microstate geometries 
that make up ‘most’ of the black hole. 

One possibility might be the generic non-BPS geometries do not have ergoregions despite the fact that they carry 
angular momentum. However, we argue that such a scenario is implausible as follows: The fuzzball description 
would now require that both the horizon and the ergosphere arise as effective surfaces in ‘coarse-graining’. However, 
quantum fluctuations must then extend out to the ergosphere. In particular, these fluctuations extend to regions 
of the spacetime which should be causally accessible to asymptotic observers on finite classical timescales. Hence 
it seems inconsistent to say that the underlying microstate geometries are hidden from asymptotic observers in this 
scenario. 

Hence as argued above, if the non-BPS microstate geometries are horizon-free with an ergoregion, they should 
expect an ergoregion instability. However, it may be that instability timescales calculated for the JMaRT solutions 
are not representative of those for typical microstate spacetimes. In particular, the latter should have complicated 
throats — as seen in their supersymmetric counterparts [13-15] — which would emulate the absorptive behavior of 
a black hole horizon. Hence it might be expected that the relevant timescales are extremely long. An important 
question is then whether the instability timescale is classically finite or not. That is, will this timescale diverge as 
the quantum numbers grow as described above. Gertainly finding more generic non-BPS microstate geometries is an 
essential step towards resolving this issue. 

In closing, we note that in the context of the AdS/GFT, a complete description has been produced for half-BPS 
microstate geometries with AdSs [55] and AdSa [56] asymptotics. This framework has given rise to an interesting 
program of semi-classical quantization [53, 57, 58] and a coarse-graining description of spacetime geometry [59]. With 
this program in mind, it is useful to recall the role of the smooth horizon-free microstate geometries in Mathur’s 
‘fuzzball’ program [11]. 

The BPS microstate geometries for the D1-D5 system can be derived by studying the Fl-P geometries and applying 
a series of duality transformations [15]. There the winding and wave numbers might be quantized by the geometry 
but classically the amplitudes of the string excitations are continuous variables. Solutions where select modes are 
excited with a large amplitude can then be seen as ‘coherent states’ of the underlying quantum theory. Such solutions 
may be further useful to understand certain properties of typical microstates, e.g., their transverse size [11]. However, 
ultimately a generic state will have a vast number of modes excited with very few quanta and hence the corresponding 
‘spacetime’ will not be accurately described by a classical geometry. However, the family of classical geometries still 
serve as a guide to the classical phase space which must be quantized [58]. 
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In the present context, we wish to go beyond the BPS sector where program is much less developed. In particular, 
we still face the challenge of constructing a more or less complete family microstate geometries. The existence of the 
JMaRT solutions indicate that at least certain non-BPS states can be described by classical geometries. However, it 
is not at all clear how large a class of nonsupersymmetric smooth horizon-free geometries exists. Going beyond the 
present special class of solutions will probably call for the development of new solution-generating techniques, but the 
JMaRT geometries offer hope that a broader class of nonsupersymmetric solutions can be found. This will certainly 
be an intriguing direction for further research and will undoubtedly lead to interesting new insights and discoveries. 
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APPENDIX A: WKB MATCHING 


In this appendix we use the usual WKB wavefunctions and WKB connection formulae at the turning points to 
relate the amplitude of the wavefunctions in the four distinct regions of the scattering problem and, in particular, to 
derive (47). The four WKB regions are (see Fig. 1): Region I, the innermost forbidden region (0 < a; < xq); Region 
II, the allowed region where V+ is below (xq < x < xi); Region III, the potential barrier region where V+ is above 
(xi < x < X 2 ); and Region IV, the external allowed region where is below V_ (a ;2 < x < 00 ). The WKB 
wavefunctions in region I and in region IV were already written in (45) and (46), respectively, and in regions II and 
III they are given by 


ffm 


Co 




exp 


i 


Vt da 




Co 




exp 


-i / y 1 ax 

J Xi 


C4 




exp 


■ TOi/> / VW\ dx 

X\ 




Cs 


my"|T|V4 


exp 




VTd. 

f V\T\dx 


(Al) 

(A2) 


Using the WKB connection formulae in each turning point, xq, xi and X2, we can find the relations between the 
amplitudes Cds (i = I, • • • , 7) of the several regions, yielding: 


C2 = CiC'^ , C3 = . 
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with 7 and rj defined in (48) and (49), respectively. Finally, combining these three sets of relations yields (47). 


APPENDIX B: DETAILS OF NUMERICAL ANALYSIS 

In this Appendix we discuss some issues related to the precision used in the numerical computations. Even though 
the very small imaginary parts of u are well described by both approximations, for completeness we show that they 
are not a numerical artifact due to loss of precision in our numeric routines or a by-product of using the approximate 
solutions to specify the boundary conditions. In Fig. 6 we plot the imaginary part of to for several values of the number 
of digits of precision used in the calculation. We use the same parameters as before and set iV = 0, I = = 4. We 
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FIG. 6: Effect of increasing digits of precision used on imaginary part of eigenvalue. 


see, as one would expect if the imaginary part were actually non-zero, that the eigenvalue converges to a constant 
value when the number of digits is larger than the size of the imaginary part. 

With only the asymptotic form of the solutions to specify the boundary conditions we are not actually setting the 
coefficient on the divergent term to zero. Instead, there will always be some amount of the divergent solution in the 
numerically defined solution. The suppression of the divergent term is dependent on how deep into the asymptotic 
region we choose to apply the boundary condition. To ensure that these small divergent terms are not causing any 
errors we study the effect of varying the point at which we apply the boundary conditions. This has been shown in 
Fig. 7 on the left and right. In both cases, we again see that the eigenvalue converges to a constant value as we 
increase the accuracy of the calculation. 
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FIG. 7: On the left we vary the point at which the inner boundary condition is applied, logj^gXo, from the value at which 
{xK?)/{(^“^/x) = 10“®. On the right, we do the same for the outer condition from the value at which {xk^)// x) = 10®. 


APPENDIX C: DETAILED ANALYSIS OF THE INSTABILITY 

The existence of a solution to the matching procedure can be reduced to the requirement that a number of constraints 
be satisfied. The difficulty one runs into when trying to discuss the general properties of these solutions is that while 
all the parameters appearing in the various equations are uniquely determined by the set {Qi, Qs, 7?, m, n}, it is 
difficult to write explicit expressions for them. In this sense, the fact that the parameters (9) can be written in such a 
simple form is really quite surprising since all are proportional to M, which can at best be defined implicitly in terms 
of the above parameters. 

Hence it is useful to have an approximation for M that allows one to understand the general behavior of the various 
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parameters. Surprisingly there is quite a simple approximate solution given by 
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where we recall that Qp = nmQiQ^/E?. For most parameter values, this expression is accurate on the order of a few 
percent. When one of the D-brane charges, say Qi, grows much larger than Q 5 ~ this approximation can break 
down, though only by a few percent times (m — (n + 1)). Similar problems appear when R? ^ Qi ^ Qs, in this 
case the error appears to be of the same order. The important thing to note is that it gives the correct scaling of M 
with the various parameters in all situations. In most cases, except those noted previously when n, it also gives 
the correct order of magnitude. Treating m as a continuous parameter, the approximation appears to produce the 
approach to the supersymmetric limit exactly. 

Using this, one can approximate, or at least bound, the parameters appearing in the solutions. 
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In the above, the inequalities result from writing sf < SiCi = QijM , in particular they become exact for the extremal 
limit. From the expression for e one sees that it is finite, and in fact positive for all m > n + 2. It is only in the 
extremal limit that e —>■ — 00 , which precludes any possible instability. One may also check from these forms that 
s/g'^ <C 1 for all values of the parameters, which can be verified numerically for sets of parameters in which the 
approximations are less trustworthy. In what follows then we will neglect e where consistent. 

The timescale of the instability is an increasing function of v which, given the above considerations, is given by 
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(C5) 


Unfortunately, we cannot make any definite statements about the size of Qp/gR^ like we did previously for e since 
it can be made arbitrarily large or small just by varying R. At this point we could use the explicit forms for Vq and 
c to discuss the general properties of the solutions. Instead we will for now set A = 0 to make the discussion more 
transparent. Non-zero A will not change the general features of the solutions. 

Setting A = 0, the above expression for n simplihes quite a bit 


(C6) 

We are now in a position to discuss the behavior of the timescale for various different solutions. Recall that the 
timescale for the instability is smallest when n is smallest. Therefore the instability will be strongest when no = I + ^ 
is smallest. This, of course, does not mean that we should necessarily consider solutions with ^ = 0, in fact we shall see 
in a moment that such solutions are not possible. More precisely the minimum value of I for which all the constraints 
can be satisfied will lead to the most unstable solution. 

Similarly, when or e/g"^ is largest the instability be the strongest. We shall deal with c next, but for now it is 
sufficient to note that it is only dependent on m and n. Observe from (C4) that for fixed R, e/g^ varies roughly 
like the inverse of the charges, therefore when one considers limits in which the charges grow, the timescale of the 
instability diverges. Similar arguments hold when R is vastly different from the charges, we find that s/g^ shrinks 
and the lowering effect of (? is diminished. It appears then that the instability will be strongest when Qi ^ Qo ^ R? ■ 
To discuss the relative effect of c we should return to the constraints. These also simplify when we set A = 0 and 
we may consider the simpler constraint c — t'o > 0. The exact form that c takes is dependent on the sign of (/. By 
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studying the constraints, it turns out that solutions with C > 0 general exist, but for larger values of I than 

when C < 0. Given the considerations above, the effect of these modes will be subdominant. We therefore focus on 

C = —nm^ + mnicj) < 0 which implies that One can then write the constraint as 

c-vq = [(to - n)(TO^ + m^) - (2A^ + 1)] - [; + 1] (C7) 

= (to — (n + l)){m,p + TO^) — {I — TO^ — m^) — 2{N + 1) > 0 . (C8) 

Further, it can be shown that when this is satisfied, the other constraints follow automatically. The last two bracketed 
terms in the final line are positive, so a solution requires that to^ + to^ > 0 , implying that must be positive. 
This is a general result that is also obtained when C > 0 or A ^ 0. When c is largest, the timescale will be shortest, 
therefore the lowest harmonic = 0 will lead to the strongest instability. One can also make c large by choosing 
and TO 0 as large as large as possible, Tto^ = I, but taking I large will not necessarily give us a very unstable mode 
because as noted before it will cause vq to rise which has an opposing effect. Since enters weighted by s/g"^, the 
more important contribution will be that from vq and the net effect is a less unstable mode. Finally, note that to^ 
and TO^ appear symmetrically in c, so that the value of v will be independent of the partition of I into and to^. 
This does not mean that the timescale will be independent of this partition since it is a weakly shrinking function of 
ICI for fixed v. When |C| is maximized the timescale will be the shortest, which is the case when = I, to^ = 0. 

To summarize then, for a fixed mode that solves the constraints, the instability will be strongest when Qi ~ ~ 

R^. On the other hand, when we fix a particular background, the instability with A = 0 will be strongest when I = 
is as small as possible and = 0 . 

Finally then we may discuss the solutions for which X ^ 0. It turns out that the various scalings of the other 
parameters appears not to be changed. When A 7 ^ 0, the constraint — Vq > 0 becomes easier to satisfy since c picks 
up a contribution proportional to dA while the contribution to t'o is smaller. The tougher constraint to satisfy is then 
the one that implies > A^. When all other parameters are fixed, this places upper and lower bounds on (i.e., we 
allow negative A) A. We will not go into detail here, but instead note one can always find solutions with non-zero A 
by going to sufficiently large angular momentum, 1. 

When studying the characteristic time for the instabilities, one finds that the timescale decreases as A is raised, but 
reaches a maximum shortly before reaching the upper bound. For negative values, on the other hand, the timescale 
is a constant decreasing function of A. As mentioned, solutions with non-zero A require larger values of I than when 
A = 0. Though larger I tends to increase the timescale, the overall effect of going to larger I to accommodate non-zero 
A can still lead to shorter timescales. 


APPENDIX D: BOUND STATES 

The general radial dependence of the scalar field at large distances from the core is determined by the sign of k^. 
When it is positive, the general solution oscillates with a power-law falloff. This is the behavior that led to the in and 
outgoing waves at infinity which we have already discussed. The other two possibilities, where is zero or negative, 
can lead to quite different behavior. For the former there is an exact solution, while the latter may again be solved 
with a matched expansion. 


1. = 0: Marginally Bound States 

By considering the special mode with = A^, both the angular and radial equations simplify sufficiently that 
an exact solution may be found. Such a choice removes all to dependence from the angular equation allowing it to 
be solved independently. The result is the eigenvalue equation for the harmonics on an S^. The exact eigenvalue is 
A = l{l + 2). 

For the radial equation, this choice of mode removes the k^x term; the same condition that previously led to the 
simplification in the near region. The previous solution in the near region (74) therefore becomes the exact solution 
in the entire spacetime. This means that asymptotically the equation has a basis of solutions in terms of 
Ignoring for now the part dependent on the KK momentum, these become r* and r~^~K These are simply the terms 
one expects from a Laplace series in four flat spatial dimensions where the angular momentum creates an effective 
radial potential. 

Asymptotic regularity requires that the component vanish whenever v > 1, leaving a field that falls off as 

The natural generalization of Friedmann’s analysis of ergoregion instabilities to five-dimensions would involve 
studying fields that fall off as r“^, therefore these modes will evade that analysis as long as > 1. The requirement 
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that removes the divergent term is similar to that for outgoing modes, except now it is an exact result 

^+IClT? = -(2iV+l) , (Dl) 

where iV is a non-negative integer. Here, however, we allow for either of the F functions in the denominator to diverge 
in eq. (76), leading to both possibilities for the sign before This is in contrast to the search for unstables modes in 
which we could neglect one of the possibilities since it was found to corresponded to ingoing damped modes. Indeed, 
since (Dl) contains terms linear in both A and w, one must consider both signs above in order to be consistent with 
the symmetry under flipping signs as in equation ( 22 ). 

In total then we have three constraints that must be satisfied for these modes. The first, = A^, fixes to to be 
an integer, meaning that there are no remaining continuous parameters characterizing the scalar field. For a general 
background then it is unlikely that the remaining constraints, in particular the one defining TV, can be solved by a 
judicious choice of the integer eigenvalues. On the other hand, fixing the set m^,TO 0 and A, there may be families of 
backgrounds for which these marginally bound states exist. 


2. < 0: Bound States 


The final possibility for solutions of the radial equation is uP < }?, or < 0. As in the case where k? is positive, 
we are unable to find an exact solution, though progress can be made through approximation. In particular, since 
the effect of the sign of n is only relevant at large distances from the core, we need only make slight modifications to 
the matched asymptotic expansion analysis presented earlier. 

To begin, we factor out the sign of by redefining k in, giving solutions that are real valued exponentials 
asymptotically. Requiring regularity therefore leaves only the exponentially damped “bound states”, localized near 
the core region. Explicitly, after having made the redefinition in (77), a convenient basis of solutions is in terms of 
modified Bessel functions of the first and second kind. 
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The first of these diverges at large x and so we require Ai = 0 for regularity. For now though, we leave Ai arbitrary, 
setting it to vanish only after we have performed the matching. 

In the matching region expanding Ii, and Ki, in the basis gives 
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Note that K^, contains both of these powers of x when expanded in the overlap region. While the contribution of the 
positive power to Ki, is relatively small, we will keep this contribution until after we perform the matching so that 
we may see how the approximate solution comes about. 

By construction, the solution in the near region (74) is unaffected by the redefinition of k. Immediately then we 
may proceed to matching the coefficients on powers of x in the overlap region. Setting A = I in the near region 
solution, we determine the coefficients Ai , A 2 in the outer region 
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As before, finding the spectrum of solutions now requires that we find values of the free parameters for which these 
equations are consistent with the boundary conditions. In particular, we now set Ai = 0 and therefore ask that 
the right hand side of (D4) vanishes. Again, rather than find such parameters numerically there is an accurate 
approximation that comes from noting that consistency requires A 2 be non-zero. This implies that the second term 
in (D4) must be non-zero and therefore any solution must come from cancellation between the two terms. Since the 
second term is suppressed by the factor k^'^, a comparable suppression must occur in the first term, again requiring the 
divergence of a T function in the denominator. This gives a quantization condition similar to that found previously 
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FIG. 8: Variation of the size of the imaginary part of w resulting from the numerical solution of (D4) as the precision is 
increased. 


Here again, the terms linear in oj and A implicit in the above equation - see definitions in eq. (17) — imply that both 
possibilities are required for consistency with the symmetry (22), though in practice both may not lead to solutions 
for which < A^. 

When V is real, this appears to give solutions for w which are purely real. Note, however, we must be careful 
in solving the constraint since, given the right combination of background charges, could become negative. For 
an arbitrary frequency in this range, Eq. (D4) will be complex so solutions where to has both real and imaginary 
parts may be possible. Such solutions cannot be found with the sort of perturbative expansion used in studying the 
outgoing modes since now it is the real part of v which gains a small correction, while the imaginary part is large. 
We can therefore no longer consider the behavior near the pole on the negative real axis defined by the real part of 
Eq. (D6). Instead, we have resorted to searching for these solutions by solving (D4) numerically. 

Generically, the root finding algorithm will produce a complex value of to that sets the equation to zero within 
a specified precision. Since the imaginary part is many orders of magnitude smaller than the real part one should 
ensure that it really is non-zero and not a numerical artifact. In Figure 8 we show the variation in the size of the 
imaginary part as a function of the tolerance used in finding the root of (D4). From this plot we see that the imaginary 
contribution is indeed just an artifact of trying to solve the complex equation. Surprisingly then it appears we can 
satisfy (D4) with a real value of w, even if that value causes v'^ < 0. That value corresponds to the solution of the 
equation resulting from taking the real part of the quantization conditions (D6). 

Since the condition (D6) is the same as for the outgoing modes, much of the analysis in Appendix C about the 
existence of solutions applies. The situation is somewhat more complicated in that one now allows modes with positive 
Lo and there are two possible solutions corresponding to the two signs in (D6), but the general characteristics of the 
solutions are the same. In particular, for the outgoing modes it was found that there are upper and lower bounds 
on the allowable values of A beyond which — A^ changes sign. In light of these bound states, we see that the full 
space of solutions may be considered as split into distinct regions based on the value of A. There is a small-|A| regime, 
in which one finds the outgoing unstable modes. This is surrounded, at larger values of |A|, by a regime where the 
bound states arise. 

This separation of the two types of modes according to the parameter A makes clear the difference in their origin. In 
particular, one can always find outgoing unstable modes that do not carry KK momentum, they need only be supplied 
with sufficient angular momentum. This is in accord with our interpretation of these solutions as the unstable modes 
predicted by Friedman which result from the existence of the ergoregion. In contrast, bound states will always result 
as long as |A| is large enough. This includes modes which carry no angular momentum, thus indicating the important 
characteristic of these solutions is their KK momentum and the effective five-dimensional mass it induces. 

Having established the existence of these bound states we should question just how close to the core region they 
are bound. The solutions are damped exponentially and so have characteristic size 
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To arrive at the final line we have used the approximation for M (Cl) found in Appendix C. The boundary of the 
ergoregion, on the other hand, is given by the vanishing of the norm of the Killing vector dt (10). We ignore the Oi, 02 
dependent contributions appearing in / to give an outer bound on the size of the ergoregion, given approximately by 
~ MCp. In terms of the variable x this means 
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nm{s ^ — s^)2cp — 


(D9) 


Whenever Qi and Q 5 are much smaller than i?^, the size of the bound state scales as Xj,/ ~ QiQb/R^ 'C 1. On 
the other hand, for large Qi and Q5 we have ~ Qi/B? where Qi is the smaller of the two. In other words, the 
size of the bound state is strongly dependent on the background. When the charges are large, the bound state will 
be mostly contained within the ergoregion, while for small charges the exponential tail of the bound state can extend 
far outside. 

Finally we can consider the possibility that the bound states have negative energy, which requires a detailed analysis 
of the energy integral (7). Examining the integrand evaluated on bound state solutions, we see that it may become 
negative near where the modulus of the scalar field peaks if the latter occurs inside the ergoregion. Though there are 
bound states for arbitrarily large values of |A|, the total energy will not be negative for all of these. Instead, the modes 
that tend to exhibit negative energy densities (in the ergoregion) only appear for a limited range of |A|, which is just 
beyond the small-|A| regime discussed above. That is, for values of A near where the ergoregion instability appears. 
When this is the case, the maximum of the modulus of the scalar field is inside the ergoregion and the phase velocity 
in the compactified direction Ey = w/A is negative, i.e., in the direction opposite to which the background is boosted. 
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